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The possible functional forms of the effective conductivity (Te of the randomly 
inhomogeneous two-phase systems at arbitrary values of concentrations are dis- 
cussed. Two different solutions for effective conductivity are found using a du- 
ality relation, a series expansion of in the inhomogeneity parameter z and 
some additional conjectures about the functional form of Ue ■ They differ from 
the effective medium approximation, satisfy all necessary requierements, and re- 
produce the known formulas for in the weakly inhomogeneous case. This 
can also signify that cTg of the two-phase randomly inhomogeneous systems may 
be a nonuniversal function, depending on some details of the structure of the 
random inhomogeneities. 



PACS: 73.61.-r, 75.70.Ak 

The electrical transport properties of the disordered systems have an impor- 
tant practical interest. For this reason they are intensively studied theoretically 
as well as experimentally. In this region there is one classical problem about the 
effective conductivity of inhomogeneous (randomly or regularly) heterophase 
system which is a mixture of N(N > 2) different phases with different conduc- 
tivities ai,i — 1,2, ...,N. We confine ourselves here by the simplest case of 
the two-dimensional heterophase systems with N — 2. Despite of its relative 
simplicity only a few general exact results have been obtained so far. Firstly, 
there is a general expression for Ce in case of weakly inhomogeneous isotropic 
medium, when the conductivity fluctuations 6a are smaller than an average 
conductivity (cr) ^ 



(1) 



where D is a dimension of the system. In our case of two-dimensional two- 
phase system (cr) = x<Ji + (1 — x)(j2, (cr^) — (cr)^ = 4a;(l — a;)(T_^, where x is 
a concentration of the first phase, cr_ = (tri — cr2)/2 , and (1) takes a form 

/ 2a_ 4a;(l - x)a- 

(Te = 0-1 1 - (1 - x) 

V CTi Zai 

a+ {l + 2{x-l/2)z-2x{l-x)z'^) , (2) 
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where C-i- = (cti + cr2)/2, and a new variable z = (7_/<t+, characterizing an 
inhomogeneity of the system, is introduced. 

The further progress in the solution of this problem is connected with the 
discovery of a dual transformation, interchanging the phases [21E1- This trans- 
formation allows to find an exact formula for CTe in the case of systems with 
equal concentrations of the phases a; = Xc = 1/2 |3] 

O-e = \/Ol(J2. (3) 

This remarkable formula is very simple and universal, since it does not depend 
on the type of the inhomogeneous structure of the two-phase system. For sys- 
tems with unequal phase concentrations a dual transformation gives a relation 
between effective conductivities at adjoint concentrations x and \ — x or in 
terms of a new variable e — x — Xc ( ^1/2 < e < 1/2 ) at e and — e 

cre(a;,cri,o-2)cre(l - x.a^.a^) = 0^02 = CTe(e, CTi, o-2)a'e(-e, cri, cr2)- (4) 

Relation (4) means that the product of the effective conductivities at adjoint 
concentrations is an invariant. Due to this relation, one can consider CTe only 
in the regions x'>Xc ( e > ) or x < a;c ( e < ) . 

However, an explicit formula for the effective conductivity at arbitrary phase 
concentrations and z attracts the main interest in this problem. One such 
formula was obtained many years ago in the so-called effective medium (EM) 
approximation ,4 , which turned out to be a good approximation for random 
resistor networks not only in the weakly inhomogeneous case . In this paper, 
using a duality relation and a series expansion in the inhomogeneity parameter z, 
we will find two explicit approximate expressions for the effective conductivity of 
two-phase systems, differing from the EM approximation. The physical models, 
corresponding to them, are introduced in other papers, where their properties 
are discussed in detail [HIIZ]. 

Let us start our investigation of the isotropic classical random two-phase 
system in the case of arbitrary concentrations with a general analysis of the 
possible functional form of the effective conductivity. Due to the linearity of 
the defining equations (mSl)i the effective conductivity of the random systems 
must be a homogeneous function of degree one of ct^ (i = 1, iV). In the case 
of = 2 , it is convenient to use, instead of Ui (i = 1,2), the variables err- 
and z (—1 < z < 1), and, instead of x, a new variable e. Then the effective 
conductivity can be represented in the following form, symmetrical relative to 
both phases, 

CTe(e,(T+,cr_) = a+/(e,cr_/CT+) = o-+/(e, z), (5) 
where ae{^^a+,a-) and f{e,z) must have the next boundary values 

a-e(l/2,cr+,cr_) = (Ti , 0-e(-l/2, (7+, cr_) = (T2, 

/(l/2,z) = l + z, /(~l/2,z) = l-z, /(e,0) = l. (5') 
The duality relation in these variables takes the form 

/(e,z)/(-e,z) = l-z2, (6) 
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from which it fohows that at critical concentration e = 

f{0,z) = VT^. (3') 

Strictly speaking, the form of a duality relation (6) is also a consequence of 
another exact relation for the effective conductivity, taking place at arbitrary 
concentrations for systems with the similar random structures of both phases 
of the system, 

(7e(e,0-i,0-2) = Cre(-e,0-2,0-l)- (7) 

This means that the effective conductivity of the random two-phase system 
must be invariant under substitution of these phases ( cri •* — > a2 ) with the 
corresponding change of their concentrations x < — > 1 — x (or e — > — e ). In the 
new variables this means that 

f{e,z) = f{-e,-z), f{-e,z) = f{e,-z). (8) 

For this reason, the duality relation can also be written in the form 

f{e,z)fie,-z) = l-z\ (9) 

It follows from (8) that the even ( ) and odd ( /„ ) parts of /(e, z) relative 
to e coincide with the even ( ) and odd ( /" ) parts of /(e, z) relative to z. 
Consequently, f{e,z) has the functional form 

f{e,z) = f{ez,e\z% (10) 

It follows from (10) that: (1) /(O, z) is an even function of z (i.e. symmetric 

in cri.2 ), (2) an expansion of /(e, z) near the point e = z = does not contain 
terms linear in e and z separately. Analogously, the odd part fa can be 
represented in the form 

Ue,z) = 2ez^{€,z), (11) 

where $ is an even function of e and z (the coefficient 2 in front of ez is 
chosen for further convenience). 

At first sight, the duality relation (11) alone is not enough for the com- 
plete determination of / in the general case. It only connects / at adjoint 
concentrations or /„ and /„ 

-fl = l- z\ (12) 

This means that /„ and fa considered at fixed z as the functions of e satisfy 
the hyperbolic relation with a constant depending on z. The relation (12) allows 
one to express /(e, z) through its even or odd parts 

/(e, z) = fa + ^/JT+T^ = fs ± Vf!-^ + z'- (13) 

For this reason, it is enough to know only one of these two parts. Usually, one 
prefers to choose an antisymmetric part as a more simple one. It follows from 
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(2) that, in the weakly inhomogeneous case, the odd part coincides with the 
odd part of (ct) and has the simplest form (compatible with (11)) 

/a(e,z) = 2e2. (14) 

As is well known, the effective conductivity in the EM approximation can be 
obtained by substitution of (14) into (13): 



aeie,z) = a+ 2ez + {2ezf + 1 - . (15) 

We will call this expression, continued on arbitrary concentrations a; = e + 1/2 
and inhomogeneties z, the EM approximation for tJe . 

However, systems with a dual symmetry usually have some additional hidden 
properties, permitting one to obtain more information about function under 
question. Moreover, in some cases these properties can help to solve problem 
exactly (see, for example, [H]). Having this in mind, we will try to investigate the 
duality relation in more detail. For every fixed z 1 (it is enough to consider 
only the region < z < 1 ) , a function / must be a monotonous function of 
e. Since the homogeneous limit z — > is a regular point of /, it will be very 
useful to expand / in a series in z 

oo 

/(e,z) = ^/fc(e)zVfc!, (16) 



where due to the boundary conditions (5') 

/o = l, /i(e)=2e. (17) 

It is worth noting here that the expansion (16) differs from the weak-disorder 
expansion of a^, in series on the averaged powers of the conductivity fluctuations 
5a/ (a) (see, for example 01). Expansion (16) is simpler, since it deals with 
variables z and e separately, while the expansion in powers of Sa/{a) is an 
expansion on the rather complicated functions of z and e. Of course, both 
expansions are connected, but expansion (16) is more convenient for our analysis 
of possible functional forms. 

Substituting the expansion (16) into (6) one obtains the following results: 

(1) in the second order on z it reproduces a universal formula (2), thus the 
latter can be considered as a consequence of the duality relation; 

(2) in higher orders there are recurrent relations between /2/c and /2/c-i, 
corresponding to connection (12); 

(3) /2fe+i(e) are odd polynomials in e of degree 2A; + 1 and /2/c(e) are even 
polynomials in e of degree 2k in agreement with (10). 

Taking into account boundary conditions (5') and an exact value (3'), one 
can show that the coefficients fk must have the next form 

/2fc+i(e) = e(l-4e2)52fc_2(e), fc > 1, 

/2fc(e) = (l-4e2)/i2fe-2(e), fc > 1, (18) 
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where g2k-2 and /i2fe-2 are some even polynomials of the corresponding degree 
and free terms of /i2fc-2 coincide with the coefficients in the expansion of (3') 

Vl - ^2 = 1 - - z^/S - zVl6 - zVl28 - zi°/256 + ... (19) 

It follows from (18) that /a is completely determined up to overall factor num- 
ber qq. Since is determined through lower fk (fc = l,2, 3) 

/4 = 4/1/3 - 3/2' = (1 - 4e2)[(85o + 12)6^ - 3], (20) 

it is also determined by the coefficient g^. Expansion (16) in the EM approxima- 
tion has a very simple form, since all g2k-2 = (A; > 1) and /2fc(e) ~ (1 — 4e^)'^ . 
Thus, we see that, in general case, the arbitrariness of / is strongly reduced 
by boundary conditions and by exact value (3) and that the third and fourth 
orders are determined only up to one constant. One can see from the EM ap- 
proximation that any additional information about function / can determine 
this constant or even the whole function. For this reason one needs to know 
what kind of functions can satisfy the duality relation (6) except general func- 
tions from (12), (13). In order to answer this question it is convenient in the 
case z 7^ 1 to pass from / to / = //\/l — z^. Then 

/(e, z)f{-e, z) = l = /(e, z)f{e, -z). (6') 

The duality relation gives some constraints on the possible functional form of 
f{e,z). For example, assuming a functional form (10), one can write out the 
next simple expression: 

/(e, z) = exp(ez(/)(e, z)), (21) 

where (/)(e, z) is some even function of its arguments. Another possible form of 
/ is 

f{e,z)^B{ez)/B{-e,z)). 

It is easy to see that they automatically satisfy eq.(6'). 

Let us now consider two simple ansatzes for a function cf). In case (a) we 
suppose that (f)(e, z) depends only on z. This means an exponential dependence 
on concentration, which sometimes takes place in disordered systems ^10) . In 
case (b) we will suppose that 0(e, z) depends only on the combination ez. This 
can signify, for example, that / depend only on mean conductivity (a) and/or 
on mean resistivity (c^^), since (cr^) ~ (l±2ez). Expanding the corresponding 
functions / in series, one can check after some algebra that it is now possible 
to determine all polynomial coefficients unambiguously! For example, one finds 
for fa in the 3-d and 5-th orders 

go = -1, 52 = -(11 + 4e^), case{a), 
go = -3, 52 = -15(1 + 126^), case{b). 
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Another way to see this is to apply boundary conditions directly to the function 
(21). In the case (a) one obtains 

^(,) = l/,lni±i, (^i±iy. (22) 

It is interesting to note that, in terms of concentration x and partial conduc- 
tivities tTj, one obtains in this case 

a, = cjlal^^. (22') 

This corresponds to the self-averaging of In cr : 

(Te = exp(cr), (cr) = x In (Ti + (1 — x) In a2, 

noted first by Dykhnc for the case of equal phase concentrations [S] and estab- 
lished later in the theory of weak localization |11J. 

In case (b) , when (j) depends only on the combination tz , one finds 

0(e.) = ^lni±^, /(e,.)^fi±^V'^ (23) 
' 2ez 1 - 2ez' ^ ' ' \l-2ez) ^ ' 

In terms of x and Ci it has the next simple form 

= VW^)- (23') 

Series expansions of (22) and (23) coincide exactly with the corresponding ex- 
pansions mentioned above. They differ from the EM approximation already in 
the third order. 

For a general form of 4){e,z), admitting a double series expansion in z^ and 

oo OO 


one can show that now /a and fi again contain one free parameter : 
go = 6((/)io — !)■ Consequently, one needs additional information or a more 
complicated ansatz for a determination of </> in the general case. This will be 
considered in another paper. 

Thus we have found two explicit functions (22) and (23), which satisfy all 
required properties. In particular, they reproduce equation (2) in the weakly 
inhomogeneous limit z <C 1. These functions can be considered as regular 
solutions of the duality relation, since they are represented by convergent series 
in z for < z < 1 except the small region z 1/2. 

The systems, having the effective conductivity just of two forms found above 
and their properties are considered in the other paper (6j (see also j7j). We give 
here only their brief description. 

The first model represents randomly inhomogeneous systems with compact 
inclusions of the second phase with finite maximal scale Im of inhomogeneities. 
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This scale can depend on concentration of the second phase Imi^ — x) (one can 
consider only the case l — x< 1/2). The stable effective conductivity ae{x,{a}) 
(here {ct} = {<7i,<72) ), depending only x and not depending on the scale on 
which the averaging is performed, can be obtained only after averaging over 
scales I > lm{x). This cre(x, {cr}) as a function of x must satisfy the next 
functional equation, generalizing duality relation (4) 

(7e(x', {a})(Te(a:" , {(t) = (jI{x, {(t}, (24) 

where x = [x' + x''^)/2. The solution of equation (24), satisfying boundary 
conditions (5'), coincides with (22) and corresponds to the finite maximal scale 
averaging approximation (FMSA) [BJ |7j . 

The second model of a random inhomogeneous systems has a hierarchical, 
two-level structure. On the first level, it consists of squares with random phase 
layers with a mean conductivity (cr) if the direction of layers is parallel to the 
applied electrical field E or with a conductivity {a^^)^^ if this direction is 
perpendicular to E. On the second level, these squares form a random parquet 
(or a lattice), which contains with equal probabilities {p = 1/2) squares with 
both orientations. Then, using universal formula (3), one can write the next 
approximate expression for 

a,(x,M) = 

which coincides with (23). 

For a comparison of the different expressions for effective conductivity (eqs. 
(15), (22) and (23)), we have constructed three plots of the corresponding func- 
tions f{e,z) at 2 = 0.8,0.95,0.999 (Fig.l) (their full 3D plots are represented 
in [7j). The lower branch in the region e > corresponds to / from (23), the 
upper branch to the EM approximation, and the middle branch to / from (22). 
It appears that all three formulas for /(e, z) , despite of their various functional 
forms, differ from each other very weakly for z < 0, 5 due to very restrictive 
boundary conditions (5') and the exact Keller-Dykhne value. This range of z 
corresponds approximately to the ratio a-ijcrx ^ 1/3. For the smaller ratios, the 
differences between these functions become distinguishable (for e > ) , growing 
significantly only for ratios 02! S 10~^. 

One can see from formulas (22), (23) that, in both cases, one gets <Te ^ in 
the limit cr2 — > , except the small region near x = 1 and z = 1. This means 
that these formulas are not valid in the percolation limit (72 — > { z —^ 1) 
for e > ^1 . One can show that such behaviour is a consequence of the 
assumptions made about the form of the function </) and/or of the structure of 
the corresponding models 

This can be connected also with a possible divergence of the series (16) in 
general case, when z — s- 1, due to a singular behaviour of cTg in the percolation 
problem piTU]. 

For this reason the formulas (22,23) cannot be applicable for the description 
of (Jef f in the limit z — > 1 (e > 0) and the corresponding percolation problem. 
It follows also from the constructed plots that EM approximation overestimates 
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Fig.l. Plots of various expressions for f{^,z) at: (a) z = 0.8, (b) z — 0.95, (c) 
z = 0.999 . 

fje [lUL I12| , and both the other formulas underestimate it in the region z 
1, e > 0. We hope to investigate this hmit in detail later. 

Thus, we have discussed possible functional forms of the effective conductiv- 
ity of random two-phase systems at arbitrary values of concentrations. It was 
shown that the duality relation and some additional assumptions about possible 
functional form of /(e, z) can give its explicit expressions, differing from EM 
approximation. They automatically satisfy the duality relation and reproduce 
all known formulas for / in the weakly inhomogeneous limit z <^ 1. 

Though the used additional assumptions are the approximate ones the ob- 
tained results (and especially an existence of the corresponding models [HIIZI) 
can be interpreted also as if ae of the two-phase randomly inhomogeneous sys- 
tems were a nonuniversal function, depending on some details of the structure 
of the random inhomogeneities. An analogous conclusion was made earlier for 
three-phase regular systems in jl3j. where a possibility to find a generalization 
of the Keller - Dykhne formula (3) for the case iV = 3 was studied numerically. 
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